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LA W OF FACILITY OF ERRORS IN TWO DIMENSIONS. 



BY E. L. BE FOREST. 
[Continued from page 9.] 

In demonstrating the important formulas (23) in the first part of this 
paper, we followed, as stated at p. 6, the analogy of a process employed by 
previous writers in finding a "multinomial formula" for the case of one va- 
riable. After the work was printed I perceived that for our present pur- 
pose it can be considerably shortened, without making it less intelligible. 

Let us dispense with notation of a, B and n as in formulas (4) and 
(5), and designate the coefficients of the first power of the polynomial by k 
as in (2), while those of the h power are represented by I. Then the given 
polynomial is 
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and its expansion to the h power is denoted by 
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Differentiate these with respect to f, and in the first of the two equations 

WSH(f). 

substitute the expressions for u, U, -=g and -^-. This gives an equation 

similar to (10), each member of which contains the product of two rectan- 
gular polynomials. Forming the coefficient of ^ i ~ 1 rf in each member, and 
equating the two to each other by the principle of indeterminate coefficients, 
we get 

If [ m '*— m,m <'(i+m),(j— "') I m *m,m '■(i— m),Q— m) j 

\ tnk—m,-iJ(i+m),(j+m) \ W^m,— m^i-m). (.;+•») 

__ (i+^)^-OT.«i ^(»+m),(j-m) ] (J m )Kn,m l(i-m),(j-m)_ (24a) 

(i-\-'m)'--m.-rJii+m),U+m) | (* m )A m ,-rJ(i->n), (j+m) 

Likewise, differentiating u and U with respect to y, and substituting in 
the second of the equations (8), 

w ® -•(£)■ ^ 

we get a second equation similar to ( 10), and equating to each other the co- 
efficients of fy' -1 in the product in each member, we have 
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The relation existing between the coefficients A and I which form the 
square groups (22), is expressed by these equations (24a) and (246). Trans- 
posing in each of them, from the second member to the first, the portion 
which does not have the coefficient i or j, and rotating each table thro' 180° 
in its own plane, and changing the signs of both members, we get the de- 
sired formulas (23). 

Now when the exponent k is made very large, or infinite, the coefficients 
I of the expansion become ordinates z to the limiting surface, and supposing 
them to be set close together so as to be consecutive, dx and dy are repre- 
sented by dx and dy, and (24) reduces to 

x = idx, y — jdy. (25) 

The square group of coefficients I in (22) contains (2m -f- l) 2 terms, and 
the whole expansion contains (2hm-\-Vf terms. Let h be made an infinity 
of the second order, as in the analogous case already cited (Analyst, Sept. 
1879, p. 147). The limiting surface will then extend over the whole of the 
infinite plane XY, and (2Am+l) 2 will represent the whole number of points 
in that plane, or the whole number of consecutive ordinates s, while (2m+l) 2 
is the number of such ordinates included in the square group under consid- 
eration; so that the portion of the plane XY occupied by this group is in- 
finitesimal, both in its length and breadth. The corresponding portion of 
the surface may therefore be regarded as an element of a plane surface, that 
is, an element of the plane which is tangent to the limiting surface at the 
middle point of the group. Let z be restricted to mean the middle one of 
this square group of ordinates, answering to l ii} in (22), and let d x z and d ;/ z 
denote the differentials of 2 taken in the x and y directions. Then the val- 
ues of all the coefficients I in (22), at the limit, may be found from z by 
successive additions and subtractions of d x z and d y z, and will stand in the 
rectangular order 

z — md x z-\-mdyZ i z+mdji-j-mdyZ ,o„> 

z — md x z — md y z \ z-fmd x z — mdyZ 
Substituting them for the corresponding values of I in (23), collecting sepa- 
rately the coefficients of z, d x z and d y z in the result, remembering that 

^— m.m * m m 1 107\ 

^j. 1 {11) 

since the sum of all the given probabilities is certainty, also observing that 
we have identically 
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and using « 1} a 2 , /9 15 /? 2 and y as auxiliary letters, thus, 
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and writing A instead of k-\-l, which is permissible because k is infinite, we 
shall find that (23) reduces to the form 

k{—o. L z-\-p x d^+r d v z ) = —i{z—a x d x z—a 2 d y z), \ ^0) 

h( — a 2 z+y dxZ + ^dyZ) = —^{z—a^z—a^dyZ). J " 

These are differential equations of the limiting surface sought, and a, /9, 
&c. are constant parameters whose significance we now inquire into. In 
all of them, the given probabilities X occur in the same order as in (2) and 
(22), and regarding them as the masses of material points in a system, and 
taking the intervals between them, Ax and Ay, as units of distance, we see 
that here as in (3), a x and « 2 are the statical moments of the system about 
the Y and X axes through the middle of the group. If therefore the given 
elementary errors are so distributed about the point of no error, which is at 
the middle of the group, as to make this point the centre of gravity of the 
system of masses A, the statical moments about any axes through this centre 
will be null, and we shall have a x — and a 2 = 0. As for ft x and /3 2 , 
the masses I in each are multiplied by the squares of their distances from the 

Y and X axes respectively, so that /3 X and /9 2 are the moments of inertia of 
the system about those , and they also represent the squares of the radii 
of gyration, since the sum of all the masses X is unity. Lastly, y is the sum 
of all the products formed by multiplying each X into its distances from the 

Y and X axes. 

It is proved in mechanics that if we have a system of masses M of mate- 
rial points in a plane, and draw any rectangular axes through their cen- 
tre of gravity, and denote by u and v the coordinates of each point M, the 
"free axes" of the system in that plane, or the principal axes through the 
centre of gravity, will be determined in position by the equation 

2I(Muv) 



ten 2<p — 



2(Mu*) — Z(Mv*)' 



(31) 
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where <p is the angle which a free axis makes with the assumed axis of the 
abscissas u. (See for instance the demonstration in Wcisbach's Mechanics.) 
If the free axis coincides with this axis of abscissas, so that f = 0, we shall 
have Z(Muv) = 0. Now y in (29) is formed in the same way as 2(Muv), 
whence it appears that if the given elementary errors are so distributed 
about the assumed X and Y axes through the point of no error, as to make 
these the free axes of the system of masses X, we shall have not only o^— 
and « 2 = 0, but y = 0. 

Of course, if the distribution of the elementary errors about the point of 
no error is entirely arbitrary, that point will not in general be the centre of 
gravity of their probabilities X, and the assumed Xand Y axes will not be 
free axes. If, however, we suppose these X and Y axes to be rotated about 
the origin, so as to become for instance parallel to the free axes, this will 
make no difference in the expansion. When the positions of the coefficients 
k in a rectangular polynomial such as (2) or (22), and those of the coeffi- 
cients I in its expansion, are referred to the same rectangular axes, these 
positions are independent of the directions of the axes so long as the origin 
remains the same. For example, let x t y t and x 2 y 2 be co-ordinates of two 
given points of error, referred to horizontal and vertical axes, and let x\y\ 
and x' 2 y' 2 be co-ordinates of the same points referred to rectangular axes 
having the same origin but making an angle with the former. Multiply- 
ing together the two corresponding terms of the polynomial, omitting their 
coefficients k, we have in the two systems 

The exponents in the second members are the coordinates of the place of 
the coefficient of the product, in the two systems. The known relations 
between such systems of coordinates as are here supposed, are 

x x =x\cos,d — y^sin 6, x 2 =a; 2 cos0 — y' 2 $ia,0, 1 ,„„, 

y t = x\s\x\ d+y^cosO, y 2 = a;' 2 sm 0-f-y' 2 cos 0, J 
and addition gives 

*i+*> = K+a' 2 )cos0— (yi+y 2 )sin0, 1 
3/i +2/2 = K +*' 2 )sin 0+{y\ +y' 2 )cos 6. I ( } 

These equations, being of the same form as (33), show that just as x^y x 
and x\y\ are coordinates of the same point in the two systems, so the point 
whose coordinates are x x -f- x 2 and y x -f- y 2 m *^ e nrst system, is the same 
point whose coordinates in the other system are x\ -f- x\ and y\ -f- y\. In 
other words, the product occupies the same position under either system. 
And this being true for all the partial products which make up the total, it 
appears that whatsoever the directions of the assumed rectangular axes may 
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be, the origin being unchanged, the places of the coefficients in the expanded 
polynomial will be the same, and consequently the limiting surface will be 
unchanged. 

Supposing now that the assumed X and Y axes have been rotated so as 
to be parallel to the free axes, we may suppose further, that while the axes 
retain this new direction, the origin is transferred to some new point, for 
instance, to the centre of gravity of the masses X. This also will make no 
difference in the positions relatively to each other, of the coefficients in the 
expansion to the k power. The only difference will be, that relatively to 
the origin, the whole system of coefficients in this expansion will be trans- 
ferred through distances, in the x and y directions, k times as great as the 
distances through which the origin was transferred. 

The change of origin amounts simply to increasing or diminishing all the 
exponents of ? in (2) by one constant number, and the exponents of rj by 
another, and this, as we have already noticed, does not alter the values of 
the coefficients in the expansion. If the expansions in the two systems are 
afterwards brought together so as to coincide, then their origins, or the 
places of exponent zero in each, will be separated by intervals, in the x and 
y directions, k times as great as the intervals between the origins in the first 
power. The case is similar to that of a polynomial of one variable, as treat- 
ed in my former article (Analyst, Sept. 1879, p. 146). It appears there- 
fore, that no matter what point is taken as an origin, or what direction is giv- 
en to the rectangular axes, the coefficients in the expansion of the given 
polynomial are always the same and in the same relative positions, and 
consequently their limiting surfaces are the same, differing only in position 
relatively to the axes. In any given case, we can choose the axes which 
are most convenient, and we shall choose the free axes of the system of 
masses X, the origin being at the centre of gravity. 

To show that our previous demonstration is applicable here, we observe 
that the units of measure Jx and Ay, which form the intervals between the 
positions of successive coefficients X in the given polynomial (2), and corres- 
pond to differences of unity in the exponents of $ and rj, may be made as 
small as we please, or even infinitesimals, without affecting the positions of 
the coefficients, provided all the exponents up to m are made large in pro- 
portion. Fractional exponents can be converted into whole numbers by 
multiplying all the exponents of f by one sufficiently large number, infinite 
if necessary, and all those of rj by another, and this will not alter the values 
of the coefficients iu the expansion. 

Our system of notation, by the way, represents the coefficient of any term 
in the first power of the polynomial by X with sub-indices equal to the ex- 
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ponents which c and rj may receive in that term. If now the 
p]ace of f <hf, is transferred from one term to another, by adding or subtract- 
ing a constant number in all the exponents of £ and another in those of tj, 
this new origin can be made the middle of a square polynomial, by adding 
to the given terms a sufficient number of terms with zero coefficients, to 
complete the square. These terms are imagined merely to show that the 
conditions of the demonstration are satisfied, and their coeffiicients A — 
disappear in the result (23). 

Suppose for instance that the poly- 1 
nomial is such that its coefficients J 
occupy the nine angular points of the I 
four small squares into which the 
square A B CD is divided in the 
marginal figure, and that we wish to 
change the positions of the axes so that I 
instead of being parallel to the sides of I 
this square, with their origin at its cen- 
tre, they shall have new directions and I 
a new origin. Describe the square E- J 
FOH with its centre at the desiredl 
origin and its sides parallel to the new axes, and large enough to include 
within it the square ABCD. We may conceive that every point of the 
plane within this new square is occupied by a coefficient X, only all of them 
are equal to zero except those at the nine angular points in the original sq. 
ABCD. All that our demonstration requires is that the coefficients X in 
the given polynomial should form a square group with one term in the cen- 
tre, the sides of the square being parallel to the axes of X and Y; and these 
conditions are fulfilled in the square EFGH. The number of terms which 
the square contains is (2m + If, and m here is of course infinite. But 
the relation (23) between these terms X and the terms I in any square group 
of equal size in the expansion will still hold true, when the products XI are 
formed in the manner there stated. The first members of the two last eq's 
in (23) will represent the moments of the system of products XI, about axes 
passing through the centre of the square group and parallel to the new axes 
of Zand Y. The numbers i and j will represent the coordinates of the 
centre of the square group of expanded terms I in the same X and F sys- 
tem, the centre of the whole expansion being at the origin, so that i and j 
may be any whole numbers, up to infinity, and idx and jdy are the coordi- 
nates themselves. But a distance, or a lever arm, is not altered by chang- 
ing the unit of measure. It makes no difference, for example, whether we 
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say 3 feet, or 36 inches, or (36Xoo) times (1 inch -f- oo). We may restore 
then, if we please, in (22) and (23), the finite units of measure Ax and Ay, 
regarding m, i and / as finite numbers, only they will now, in general, be 
fractional. The same is true of the m in (26), for d x z and d y z are the incre- 
ments of z corresponding to those increments in x and y respectively which 
are equal to the units of measure Ax and Ay, and which in the limiting sur- 
face are reprseented by dx and dy. Thus the validity of the differential eq's 
(30) of the limiting surface remains unimpaired by the change of axes, 
provided it is understood that the moments a, /3, y of the given system of 
masses X, as expressed in (29), are now to be taken with reference to the new 
axes, that is, the free axes of the system, while the coordinates x and y, or 
idx and jdy, refer also to those axes. 

Adopting therefore the centre of gravity of the coefficients X as the origin, 
and their free axes as the axes of X and Y, we have in (29), as before sta- 
ted, fltj, a 2 and y equal to zero, and (30) takes the simple form 

kft x d x z = — iz, kft 2 d y z = — jz. (35) 

Giving to i and j their values from (25) we can obtain 

4i! = — xdx dji __ — ydy ,„„, 

z kptidx) 2 ' z kP»{dy)*" { ' 

According to what was said in connection with (29), ^ x (dxf and fl^dy) 2 
represent the squares of the radii of gyration of the masses X about the X 
and Faxes, when Ax and Ay are represented by dx and dy at the limit. 
Let us write 

r\ = p x {dxf, r\ = p 2 {dyf. (37) 

By the theorem which I gave in Analyst, May 1 880, p. 78, lcr\ and 
hr\ represent the squares of the radii of gyration of the coefficients I in the 
whole expansion to the h power, so that r x y/h and r 2 j/7c are the "quadratic 
mean errors" in the x and y directions. It is immaterial to the validity of 
that theorem, whether the radii of gyration are taken with reference to axes 
passing through the centre of gravity and parallel to the X and Y axes 
originally assumed, or are taken with reference to the free axes which stand 
at an angle <p with the others, for if the positions of the coefficients in the 
first power of the polynomial and its expansion are referred to any rectan- 
gular axes with origin at the centre of gravity, it will always hold true, as 
in my article cited, that the squares of the radii of gyration in the h power 
are h times the the squares of the corresponding radii in the first power. 
In that demonstration the exponents m, n, p and q may be as large as we 
please, even to infinity, and the units of measure Ax and Ay as small as we 
please, or infinitesimals. Whatever the direction of the axes may be, the 
polynomial can always be put in the form of a rectangle with sides parallel 



—48— 

to the axes, by adding a sufficient number of terms having zero coefficients, 
as we have already noticed. 

It is a known mechanical property of the free axes, that the radius of 
gyration is a maximum about one free axis and a minimum about the other. 
If then a given polynomial is raised to successive powers, making k = 2, 
& = 3, & = 4, &c, and the coefficients in these expansions are all referred 
to the same axes passing through the centre of gravity of the coefficients in 
the first power, the radii of gyration in the k power are \/k times what they 
are in the first power, for any given position of the axes, so that they are 
still a maximum when taken about one of the free axes of the first power, 
and a minimum abeut the other. Hence, the free axes retain a constant 
position, in all the successive powers. 



If we now write 



A ^2^f' ** = **¥ W 



d J? = — 2h\xdx, ^ =—2hlydy. (39) 

z z 



the equations (36) become 

dji 
z 

Their integral is 

z = ce- ( * lte8+)W >. (40) 

They are readily derived from it by differentiating it with respect to x and 
y and dividing by z. To find the value of c, we consider that as the sum of 
all the coefficieuts X in the first power of the polynomial is unity, the sum 
of all the coefficients z in its expansion to the k power must be unity. 
Hence arises the condition 

-J J zdxdy=l, (41) 



which gives 



dxdyJ _ m - 



c 



e-^dxj e-w>dy = 1. 



dxdy* 

The known values of the two definite integrals are -[/tz-^-K^ and ^/n-~h 2 , so 
that we get 

and the complete equation of the limiting surface stands 

hxh^dxdy ^^^a+h^t) ,^\ 
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[To be continued.] 



